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The QCD partition function for the Wilson Dirac operator, Dw, at nonzero lattice spacing a can 
be expressed in terms of a chiral Lagrangian as a systematic expansion in the quark mass, the 
momentum and a 2 . Starting from this chiral Lagrangian we obtain an analytical expression for 
the spectral density of Ys{Dw +m) in the microscopic domain. It is shown that the 75-Hermiticity 
of the Dirac operator necessarily leads to a coefficient of the a 2 term that is consistent with the 
existence of an Aoki phase. The transition to the Aoki phase is explained, and the interplay of the 
index of Dry and nonzero a is discussed. We formulate a random matrix theory for the Wilson 
Dirac operator with index v (which, in the continuum limit, becomes equal to the topological 
charge of gauge field configurations). It is shown by an explicit calculation that this random 
matrix theory reproduces the a 2 -dependence of the chiral Lagrangian in the microscopic domain, 
and that the sign of the a 2 -term is directly related to the 75-Hermiticity of Dw- 
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1. Introduction 

For sufficiently small eigenvalues the spectrum of the QCD Dirac operator can be obtained 
from the zero momentum part of a chiral Lagrangian that is consistent with the global symmetries 
of QCD [1]. The reason is spontaneous symmetry breaking and the existence of a mass gap, so that 
QCD at low energy is described by a chiral Lagrangian for the Goldstone bosons. If the Goldstone 
bosons corresponding to the mass scale of the Dirac spectrum we are interested in, have a Compton 
wavelength that is much large than the size of the box, the zero momentum part of the partition 
function, which determines its mass dependence, factorizes from the nonzero momentum part [2]. 
All theories with a mass gap and the global symmetries and symmetry breaking pattern of QCD 
are described by the same zero momentum chiral Lagrangian. Analyzing the simplest such theory, 
chiral random matrix theory, allows us to answer questions that cannot be addressed otherwise. 

The Dirac operator in lattice QCD has discretization errors, and if the discretization errors 
change the underlying symmetries, the low lying Dirac spectrum will be affected significantly. 
This is the case for the Wilson Dirac operator D = Dyj + m, where 

^ = ^(V M + V*)-iaV*V M (1.1) 

is written in terms of forward (V^) and backward (V* ) covariant derivatives, and m is the quark 
mass. The Wilson term violates the axial symmetry but the Dirac operator remains 75 Hermitian 

D t w = 75Dw75- (1-2) 

Because of this, the operator D5 = Ys(Dw + m) is Hermitian, and is computationally much simpler 
to work with. Its spectrum has been analyzed in lattice QCD simulations [3, 4] and by means of 
mean field studies of chiral Lagrangians [5, 6]. 

For a = 0, the spectrum of D5 has a gap [—m,m]. At nonzero a, states intrude inside the gap, 
and for sufficiency large a, the gap closes, and we enter what is known as the Aoki phase (see Fig. 
1). In terms of Dw, this is the point where m hits the cloud of eigenvalues (see Fig. 1, right). 




Figure 1: Schematic picture of mean field results for the specttal density of D$ (left two figures) and the 
spectrum of Dw (left). The right two figures show the onset of the Aoki phase. In this Figure, W = Ws/E. 



The sharp edge of the spectrum of D5 in the mean field limit smoothes out for finite V. In lattice 
QCD simulations [3], a roughly Gaussian tail inside the gap with a width that scales as a/^fv is 
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observed. We will provide an analytical explanation of this result by means of a chiral Lagrangian 
for the Wilson Dirac spectrum. We will also construct a random matrix theory with the symme- 
tries of the Wilson Dirac operator at a ^ that reproduces the zero momentum part of the chiral 
Lagrangian for the Wilson Dirac operator including 0(a 2 ) terms. 

2. Chiral Lagrangian 

The chiral Lagrangian for Wilson fermions is an expansion in the mass, the momentum and the 
lattice spacing. Up to a few low-energy constants it is determined uniquely by the transformation 
properties of the QCD partition function [7, 8, 9]. We decompose the partition function as [10] 

Z(m,z;a) =Y i Z% !f (m,z;a). (2.1) 

V 

In the microscopic domain, where mV, zV and a 2 V are kept fixed for V — > °°, the partition function 
with index V is given by (we will see below that v is the index of the Dirac operator) 

ZX- V (m, Z ;a) = { df /de tV f / e V[i(m +Z )ZT rJ / + I(m- Z )ET rJ /t_ a 2 W8Tr(f/2+c/ t^ )] _ 
f Ju(N f ) 

Here, £ is the chiral condensate, W% is a low-energy constant and z is the axial mass. Terms that 
can be expressed as squares of traces of U and are possible [11], but they are suppressed in the 
the large A^-limit [12] and are not considered in this talk. It is instructive to perform a mean field 
analysis of this chiral Lagrangian valid in the thermodynamic limit at fixed m and a. It turns out 
that the properties of the partition function depend crucially on the sign of Wg. 
W% > 0: For z = the partition function shows a second order phase transition at mL = Sa 2 Wg. At 
this point the quark mass hits the cloud of eigenvalues of Dw- The spectrum of D5 follows from the 
z-dependence of the partition function. It has a gap [— z g ,z g ] that closes at the critical point, where 

^ = m (l-(8a 2 W 8 /mr) 2 / 3 ) 3 / 2 . (2.3) 

W& < 0: For z = the partition function does not undergo a phase transition. This implies that the 
eigenvalues of Dw are not scattered in the complex plane which is consistent with an anti-Hermitian 
Dirac operator, and violates the very premise we started form. One possibility is that Dw would be 
both anti-Hermitian and 75-Hermitian, but then {ys,Dw} = 0, which is not the case. 

2.1 QCD inequalities and the sign of Wg 

A direct consequence of D f 5 = D5 is the QCD inequality 

Z^ 2 v (m,z) = (det 2 (Y 5 (D w + m)+z))>0 for m,z real. (2.4) 
For an anti-Hermitian Dirac operator we have that 

r 2d ^Z^(m, z )=r 2d]m ^(det 2 ({D w + im) + iz))>0 for m,z real, (2.5) 
where dim(D) is the total dimension of the Dirac matrix. In the continuum limit, / 2dim ( D ) = (— l) v . 
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By changing variables U — > iU in the chiral partition function (2.2) it follows that 

Z*^(0,0,W 8 ) = (i) N f v Z*' Nf (0,0-Wz). 

Since for large mass the partition functions for +Wg and — W$ have the same sign, one of them must 
change sign as a function of m, and not both signs of Ws can be allowed by the QCD inequalities. 
In Fig. 2 we show the Nf = 2 partition function for v = 1 as a function of m (left) and as a function 
of im (right). The blue curves show the W% > results, whereas W% < for the red curves. 
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Figure 2: The mass dependence of the partition function (2.2) for real mass and imaginary mass. 

The conclusion is that W% > can only correspond to a 75-Hermitian Dirac operator, whereas 
for W% < the chiral partition function only has a definite sign for imaginary mass so that the 
corresponding Dirac operator must be anti-Hermitian. 

3. Dirac Spectrum 

To find the Dirac spectrum of D5 at z we have to add a fermionic quark with axial mass z and 
a bosonic quark with mass z' to the partition function so that the generating function is given by 



Z QCD, 



(m,m ,z,z \a) 



I det( / 5(% + m)+z) -TTdeKzy +m f )\. 

\ det( 75 (D w + m')+Z / -/£) A f A 1 W fJ 



(3.1) 



The infinitesimal ie term is essential. The resolvent 



dz! 



= ( Tr— l - V\det( D w + m f ) 

z < =z \ Y5(Dw+m)+z-ie L / y w *' 



(3.2) 



is discontinuous across the support of the spectrum, and for real z its imaginary part gives the 
spectral density of the Dirac operator (denoted by p^(z)). The chiral Lagrangian corresponding 
to the generating function (3.1) is uniquely determined by symmetries. What is not determined 
by symmetries is the integration contour of the group manifold. For compact integrals this is not 
an issue, but for noncompact integrals corresponding to bosonic quarks, a change in integration 
contour may lead to a different analytical continuation of the integral in its complex parameter 
plane [13, 14, 6]. When Dw is Hermitian or anti-Hermitian it is possible to analytically continue 
the partition function in the entire complex m-plane. When Dw is non-Hermitian the generating 
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function is not analytic if m is inside the support of the spectrum and the expression (3.2) cannot 
be used to derive the Dirac spectrum. 

For a 75-Hermitian Dirac operator the partition function (3.1) is an analytical function of z, and 
only the part of the complex a-plane to which the partition function can be continued analytically 
can correspond to lattice QCD at a 7^ 0. In particular, analyticity in z requires that Re(a 2 W^) > 0. 

The second interesting observable that can be calculated from the generating function is 



£ v (m) = logZ v (m,m ,z,z;a) 
am 



(3.3) 

z'=z=0,m'=m 



For small a it can be shown that the discontinuity of this quantity accross the real axis gives the 
distribution of the chirality over the real eigenvalues of Dw, 

p v x (X w ) = l £ S(A w -Af)sign(<% 5 |£>)Y (3.4) 

WfeR / 

where are the eigenvalues of Dw- The integral J dX w p^(X w ) is equal to the index of the Dirac 
operator. It is defined configuration by configuration as £ t ^w €K sign((A:|75|A:)). 

In the microscopic domain of QCD, the generating function (3.2) is given by an integral over 
the zero momentum part of the Goldstone manifold 

Z^ SY ' v (m,m',z,z';d) = jdU sdet(^) v e^ STr ^[ c/ - c/ "]) + ^ STr (^[ f/+c/ "» + ' j2sTr ( c/2+c/ ").(3.5) 

The integration is over the maximum Riemannian graded submanifold of G/(l|l) [13] and can be 
worked out analytically. The result is a sum of products of simple one-dimensional integrals [10]. 
Calculating the integral / dX w p^(X w ) we find that it is identically equal to v. This shows that the 
parameter v in the chiral Lagrangian is the index of the Dirac operator defined below Eq. (3.4). 

It is straightforward to extend this result for QCD with dynamical quarks. The result for Nf=l 
is discussed in [15] in this volume. 

4. Random Matrix Theory 

Since the chiral Lagrangians (2.2) and (3.5) are determined uniquely by the global symmetries 
of QCD, they can also be obtained from a random matrix theory with the same symmetries. The 
random matrix partition with index v is defined by [10] 



yRMT. v 



= J dAdBdW det N f(D w + m + zY5) P(D W ), (4.1) 



with 

Here, A is a square matrix of size nxn, and B is a square matrix of size in + v) x [n + v). The 
matrix C is a complex n x {n + v) matrix. The matrix elements of Dw are taken to be distributed 
with Gaussian weight 

P(A,B,W) = e -^+#]-Pr[wwn^ (42) 
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where N = 2n + v. For a — >• the random matrix Dirac operator Dw has exactly v zero eigenvalues. 
For nonzero a, the index of the Dirac operator defined below Eq. (3.4) is exactly equal to v for 
each matrix in the ensemble. 

In the microscopic domain, this random matrix theory partition function reduces to the chiral 
Lagrangian introduced in Eq. (2.2) with L = N/V and Wg = \a 2 N. Notice that Wg is positive. 



5. Results for the Dirac Spectrum 



Analytical expressions for pJ(As) and p^(Xw) were derived from (3.5) in [10]. Plots of p^X^) 
are shown in Figs. 3 and 4. For a = the eigenvalue density of D5 can be decomposed as 



p 5 v (A 5 ) = vS(X 5 -m) + p Xs>m (X 5 ). 



(5.1) 



For a ytz the width of the peak at x = m becomes finite. In Fig. 3 we show the spectral density of 
D5 for v = 1 and a = 0.05. Also shown is the a 2 W%V <C 1 analytical result for the density of the 
"topological" eigenvalues for V = 1 which is given by a simple Gaussian 



VE 2 (A 5 -m) 2 



(5.2) 



For small a, it can be shown for arbitrary V that this distribution approaches the spectral density of 
the real eigenvalues of Dw up to a shift by m [11]. This figure suggests that it is possible to extract 
the low-energy constant Wg from the width of the distribution of the smallest eigenvalue. 













rh = 3, a = 0.05 
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Figure 3: Microscopic spectral density of 
75 (Dw +m) for v = 1 (blue curve). Also 
shown is the distribution of the smallest eigen- 
value for v = 1 obtained in the a — > limit. 



Figure 4: The microscopic spectrum of 
75 (Dw+m) formVE = 3, V = 2 and a^/W%V = 
0.125, 0.250 and 0.500, respectively. 



For | A5 — m 1 1 a fixed for a — > the spectral density inside the gap can be obtained from a saddle 
point analysis: 

p 5 (X 5 )~e- z2v ^- m f/ 16alw * for 0<x«m. 

The width parameter is given by a 2 = 8a 2 Wg/V£ 2 , and the ratio of a and and AA = TlfLV is 
equal to a /AX = /n)ay/W&V [10]. This is exactly the scaling behavior found in lattice QCD 
simulations [3]. We note that in this reference the results are an average over different values of 
V. Our results show that, in particular for small a, the Wilson Dirac spectrum shows a strong 
dependence on v, and it would be very interesting to analyze the lattice results for fixed index of 
the Dirac operator. 
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6. Conclusions 

Using chiral perturbation theory we have obtained exact analytical expressions for the spec- 
trum of the Wilson Dirac operator in the microscopic domain. Our result are obtained for fixed 
value of the index of the Dirac operator and corrections up to 0(a 2 ) have been included. We have 
formulated a random matrix theory that in the limit of large matrices reproduces the chiral La- 
grangian and the analytical results for the Dirac spectrum. Using QCD inequalities for the Wilson 
Dirac operator we find that the sign of the coefficient of 0(a 2 ) in the chiral Lagrangian is consistent 
with the existence of an Aoki phase. In the limit of small a we have obtained analytical results for 
the distribution of the smallest Dirac eigenvalues with fixed nonzero index of the Dirac operator. 
We find that the width of the distribution scales as a/y/V in agreement with earlier lattice simu- 
lations. Our results make it possible to extract additional low-energy constants of Wilson chiral 
perturbation theory from the distribution of the smallest Dirac eigenvalues. 
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